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Abstract 

Investigating thermodynamic properties of liquid-solid transitions of flex- 
ible homopolymers with elastic bonds by means of multicanonical Monte 
Carlo simulations, we find crystalline conformations that resemble ground- 
state structures of Lennard- Jones clusters. This allows us to set up a 
structural classification scheme for finite-length flexible polymers and their 
freezing mechanism in analogy to atomic cluster formation. Crystals of 
polymers with "magic length" turn out to be perfectly icosahedral. 
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1 Introduction 

Small crystals of cold atoms such as argon [1] and spherical virus hulls enclosing 
the coaxially wound genetic material [3J[3] exhibit it, and - as we will show here - 
also elastic flexible polymers in the solid state: an icosahedral or icosahedral-like 
shape. But why is just the icosahedral assembly naturally favored? The reason 
is that the arrangement of a finite number of constituents (atoms, proteins, 
monomers) on the facets of an icosahedron optimizes the interior space filling 
and thus reduces energy. 

Many-particle systems governed by van der Waals forces are typically de- 
scribed by Lennard- Jones (LJ) pair potentials 

E w (r ij ) = 4e[(a/r ij ) 12 -(a/r ij f}, 
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Figure 1: (a) Anti-Mackay and (b) Mackay growth over layer on the facet of an 
icosahedron [5]. 

where is the distance between two atoms located at and (i, j — 1, . . . , N), 
respectively. An important example are atomic clusters whose structural prop- 
erties have been subject of numerous studies, mainly focusing on the iden- 
tification of ground-states and their classification. It has been estimated (4j 
that icosahedral-like LJ clusters are favored for systems with N < 1690 atoms. 
Larger systems prefer decahedral structures until for N > 213000 face-centered 
cubic (fee) crystals dominate. In the small-cluster regime, energetically optimal 
icosahedral-like structures form by atomic assembly in overlayers on facets of 
an icosahedral core. There are two generic scenarios (see Fig.[l|: Either hexag- 
onal closest packing (hep) is energetically preferred (anti-Mackay growth), or 
the atoms in the surface layer continue the fec-shaped tetrahedral segment of 
the interior icosahedron (Mackay growth) [5 J. The surprisingly strong depen- 
dence of structural liquid-solid transitions on the system size has its origin in the 
different structure optimization strategies. "Magic" system sizes allow for the 
formation of most stable complete icosahedra (N = 13, 55, 147, 309, 561, 923) [BJ. 
Except for a few exceptional cases — for 13 < N < 147 these are N — 38 (fee 
truncated octahedron), 75-77 (Marks decahedra), 98 (Leary tetrahedron), and 
102-104 (Marks decahedra) [4] — LJ clusters typically possess an icosahedral core 
and overlayers are of Mackay (N = 31-54,82-84,86-97,99-101,105-146,...) or 
anti-Mackay (N = 14-30, 56-81, 85, . . .) type [7]. 

Although it seems that there are strong analogies in liquid-solid transitions 
of LJ clusters and classes of flexible polymers, comparatively few systematic 
attempts were undertaken to relate for finite systems structural properties of 
LJ clusters and frozen polymers 019]. In contrast, the analogy of the generic 
phase diagram for colloids and polymers has been addressed in numerous studies 
(for a recent overview, see, e.g., Ref. [TU]). 

Here, we systematically analyze for a frequently used flexible LJ polymer 
model with anharmonic springs [111 I12j the formation of frozen conformations in 
conformational liquid-solid transitions and find that for the chain lengths studied 
(N < 309) the ground-state conformations resemble structures known from LJ 
clusters of corresponding size. Our results for the peak structures of energetic 
and structural fluctuating quantities and the interpretation of the liquid-solid 
transitions invalidate recent results from studies of the same model [13], but 
are consistent with former studies of homopolymers with harmonic bonds [9] 
and analyses of LJ cluster formations [7 ]. In analogy to a recent analysis of 
the crystallization of flexible lattice polymers [15] and LJ cluster studies [7], 
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there is no clear size-dependent scaling behavior for energetic and structural 
fluctuations near the liquid-solid transition in the icosahedral regime being in 
the focus of our study. As expected, liquid-solid and 6 collapse transition seem 
to remain separate transitions in the thermodynamic limit [101 1151 116j . 

2 Model and method 

For our study, we employ a polymer model with truncated and shifted pairwise 
LJ potential, 

KTVtf) = E hJ (mm( rij ,r c )) - E u (r c ), 

and finitely extensible nonlinear elastic (FENE) anharmonic bonds [IT] [12] be- 
tween adjacent monomers, 

E FENE (ni+i) = -KR 2 ln{l - [(r ii+1 - r Q )/R] 2 } 1/2 - 

For the parametrization we follow Ref. |12) . In our units, the LJ parameters 
are set to e = 1 and a = 2~ 1 / 6 ro with the potential minimum at ro = 0.7 and 
the cutoff at r c = 2.5c The FENE potential has a minimum coinciding with 
ro and it diverges for r — > tq ± R with R = 0.3. The spring constant K is set 
to 40. Finally, the total energy of a polymer conformation X = (ri, . . . , r^v) is 
given by 

N N-l 

»,j=i i=i 

In order to sample the entire state space with high accuracy, we performed 
multicanonical Monte Carlo simulations [17j[T8], where recursively calculated 
multicanonical weight factors [HJJ [20] deform the energy histogram in such a 
way that entropically suppressed low-energy conformations are sampled as fre- 
quent as high-energy states which are of less interest for the conformational 
transitions discussed here. In the final production runs, at least 10 10 updates 
were generated, 90% of which were simple displacements of single monomers and 
10% bond exchange updates altering monomer linkages. All conformations we 
refer to as ground states in the following were identified in these multicanonical 
simulations and refined with standard minimizers. 

3 Results and discussion 

Before discussing size-dependent thermodynamic properties of structure forma- 
tion in liquid-solid transitions of the flexible FENE polymers, we discuss the 
structural composition of ground-state conformations in the icosahedral regime. 
Since the FENE bond potential enables almost nonenergetic deformations of 
covalent bonds near the LJ minimum, polymer ground-state morphologies re- 
semble corresponding shapes of LJ clusters. It should be noted, however, that 
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Figure 2: "Magic", i.e., perfectly icosahedral ground-state conformations of 
elastic flexible polymers with 13, 55, 147, and 309 monomers, as found in the 
multicanonical simulations. The characteristic Mackay shells are shaded differ- 
ently. 

the covalent bonds explicitly break rotational symmetries and, in consequence, 
the polymer ground states are highly degenerate and metastable. 

At zero temperature, similar to LJ clusters, also flexible FENE polymers 
with N > 13 monomers usually build up structures with icosahedral cores. The 
smallest icosahedral structure can be formed by the 13mer: A central monomer 
is surrounded by twelve nearest neighbors. For larger system sizes, additional 
monomers form new shells around the central 13mer and the next complete 
icosahedron with N = 55 contains two completely filled shells — and represents 
the core of the triple-layer icosahedron with N — 147 etc. Corresponding 
ground-state conformations with almost perfectly icosahedral shape as found 
in our multicanonical simulations are shown in Fig. [5J In complete analogy 
to the LJ clusters, ground-state morphologies of FENE polymers with inter- 
mediate chain lengths "grow" by forming anti-Mackay or Mackay overlayers, 
except for special cases such as the polymer with N = 38. If only a few excess 
monomers remain after compact core formation, these are bound to the facets 
of the icosahedral core in the optimal distance, thus filling the available space 
at the facet most efficiently by forming an anti-Mackay layer [Fig. QJa)]. For 
larger systems, a crossover to the formation of Mackay-type layers [Fig. [ljb)] 
occurs as it allows to place additional monomers on the edges. This increase in 
structural compactness energetically overcompensates the occurrence of partly 
nonoptimal distances which induces strain into the structure (this is the reason 
why for very large systems fee crystals are favored [4]). 

The liquid-solid transition of the rather short polymers considered here is 
not a phase transition in a strict thermodynamic sense and its strength depends 
on the chain length. Figure [3Ja) shows the specific heat Cy = {(E 2 ) — (E) 2 )/T 2 
as a function of the temperature T for several exemplified polymers with dif- 
ferent chain lengths, and in Fig. (3Jb), the fluctuations d(r gyl )/dT of the radius 
of gyration r gyr = \/J2i( r i ~~ r) 2 /N (with f = J2i r i/N) are plotted. The lat- 
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Figure 3: (a) Specific heats (in the inset shifted by a constant value for the 
sake of clarity) and (b) fluctuations of the mean radius of gyration d(r gyr ) /dT 
as functions of the temperature for chains of different length. Jackknife [HI [22] 
error bars are as small as the line thickness. 



ter curves exhibit noticeable fluctuations in the high-temperature regime which 
indicate the collapse transition from extended random coils to more compact, 
globular conformations (the "liquid" state). Energetic fluctuations are rather 
weak such that this transition can hardly be identified in the specific-heat curves. 
Much more striking is the sharpness of the low-temperature peaks exhibiting 
clear signals for a further structural compactification [Fig. [3]Jb)] and energetic 
optimization [Fig. 02[a)] — the liquid-solid transition^- 

There are three characteristic features: First, the transitions are particu- 
larly strong for chains with "magic" length, which possess perfectly icosahedral 

1 Our results are in correspondence with studies of L J clusters [7] and polymers with har- 
monic bonds [9]. However, there is significant mismatch with the former FENE polymer melt- 
ing studies of Ref. |13M14j . where the energy space was artificially restricted to E £ [—AN, 0), 
thus neglecting the lowest-energy states, which are, however, statistically relevant for the 
thermodynamic analysis of the liquid-solid transition. 
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Figure 4: Smallest single-shell cores of compact conformations: (a) icosahedron, 
(b) elongated pentagonal pyramid, (c) cuboctahedron (fee), and (d) incomplete 
icosahedron. Sticks represent nearest-neighbor contacts in the shell, not bonds. 

ground-state morphology (e.g., N = 13,55,147,309). A second type of liquid- 
solid transition consists of two steps, at higher temperatures the formation of 
an icosahedral core with anti-Mackay overlayer that transforms at lower tem- 
peratures by monomer rearrangement at the surface into an energetically more 
favored Mackay layer ( "solid-solid" transition) . This is the preferred scenario for 
most of the chains with lengths in the intervals 31<7V<54or81<7V<146 
that make the occupation of edge positions in the outer shell unavoidable. In 
most of the remaining cases, typically anti-Mackay layers form. To conclude, 
although the elastic polymers are entropically restricted by the covalent bonds, 
we find that the general, qualitative behavior in the freezing regime exhibits 
noticeable similarities compared to LJ cluster formation. 

For the quantitative analysis of the structuring process we introduce an 
"order parameter" enabling the discrimination of conformational phases in the 
different scenarios. At low temperatures, all monomers in the core are always 
surrounded by 12 almost equidistant nearest neighbors. However, the arrange- 
ment of these neighbors can differ. The lowest-energy assembly is the smallest 
icosahedral cell with 30 nearest- neighbor contacts in the shell [Fig. HJa)]. All 
icosahedral conformations have a central icosahedral cell like this, but it is also 
found in noncentric parts of icosahedral clusters with a sufficiently large anti- 
Mackay overlayer. The elongated pentagonal pyramid shown in Fig. |4jb) pos- 
sesses only 25 contacts in the shell and is found in structures with icosahedral or 
decahedral symmetry. The cuboctahedron depicted in Fig. HJc) is an example 
for a cell with 24 shell contacts and is the seed of compact fee geometries but 
also occurs in larger structures, since any LJ icosahedron can be considered as a 
combination of 20 fec-tetrahedra [BJ. Near the surface, also incomplete icosahe- 
dra with 11 nearest neighbors forming 25 contacts are found [Fig. [4jd)] , which 
we will not distinguish from complete icosahedral cells in the following. 

Based on these features, icosahedral cells in a given conformation are easily 
identified in its contact map. We consider two monomers i, j as being in contact 
if rij < 0.8 which is slightly larger than the optimal LJ distance r . Eventually, 
the total number of icosahedral cells n- lc is an appropriate parameter to deter- 
mine the overall geometry, at least at low temperatures: If n- lc > 1, the cluster 
is in an icosahedral state with a larger anti-Mackay overlayer, while for n; c = 1 
the shell is of Mackay type. For rii c = 0, the cluster is non- icosahedral (e.g., 
fee- like or decahedral). 
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Figure 5: Temperature dependence of the probability of icosahedral (Mackay: 
n- lc = 1, anti-Mackay: n\ c > 1) and nonicosahedral (ri; c = 0) structures for 
exemplified polymers with lengths N = 31,38, and 55. In the left panel, the 
corresponding lowest-energy morphologies are shown. 

In the following, thermodynamic properties of structure formation are ana- 
lyzed for exemplified polymers (N = 31, 38, 55) with respect to the icosahedral 
content. These examples are representatives for the different relevant structur- 
ing scenarios accompanying the crystallization. Figure [5] shows the respective 
populations p for the three structural morphologies, parametrized by Ti[ c , as a 
function of temperature. 

For N — 31, we find that liquid structures with rii C — dominate above 
T = 0.5, i.e., no icosahedral cells are present. Decreasing the temperature 
and passing T ks 0.4, nucleation begins and the populations of structures with 
icosahedral cells (n ic = 1 and n lc > 1) increase. The associated energetic 
fluctuations are also signaled in the specific heat [see inset of Fig. [3] (a)] . Cooling- 
further, n lc is always larger than 1, showing that the remaining monomers build 
up an anti-Mackay overlayer and create additional icosahedra. Very close to the 
temperature of the "solid-solid" transition near T w 0.04, where also the specific 
heat exhibits a significant peak, the transition from anti-Mackay to Mackay 
overlayer occurs and the ensemble is dominated by frozen structures containing 
a single icosahedral cell surrounded by an incomplete Mackay overlayer. 

The exceptional case of the 38mer shows a significantly different behavior. 
A single icosahedral core forms and in the interval 0.08 < T < 0.19 icosahedra 
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with Mackay overlayer are dominant. Although the energetic fluctuations are 
weak, near T w 0.08, a surprisingly strong structural crossover to nonicosa- 
hedral structures occurs: the formation of a maximally compact fee truncated 
octahedron. 

Finally, the "magic" 55mer exhibits a very pronounced transition from un- 
structured globules to icosahedral conformations with complete Mackay over- 
layer at a comparatively high temperature (T « 0.33). Below this temperature, 
the ground-state structure has already formed and is sufficiently robust to resist 
the thermal fluctuations. 

4 Concluding remarks 

In this study, we have precisely investigated thermodynamic properties of the 
liquid-solid transition of elastic flexible polymers of finite length and demon- 
strated that the formation of icosahedral and nonicosahedral structures pro- 
ceeds along similar lines as for clusters of atoms governed by van der Waals 
forces. Morphologies of polymer ground-state conformations coincide with cor- 
responding LJ clusters. However, covalent bonds between adjacent monomers 
reduce at higher temperatures the entropic freedom and thus the behavior of 
polymers in both the liquid and random-coil phase is different. Introducing a 
structural "order" parameter probing the icosahedral content of a structure, we 
have also analyzed the nucleation process in detail, which typically starts with 
the formation of an icosahedral seed, even for polymers preferring nonicosahe- 
dral lowest-energy conformations. 
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